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ABSTRACT 

We present the color distributions of globular cluster (GC) systems for 100 Virgo cluster early- 
type galaxies observed in the ACS Virgo Cluster Survey, the deepest and most homogeneous survey 
of this kind to date. While the color distributions of individual GC systems can show significant 
variations from one another, their general properties are consistent with continuous trends across 
galaxy luminosity, color, and stellar mass. On average, galaxies at all luminosities in our study 
(—22 < Mb < —15) appear to have bimodal or asymmetric GC color distributions. Almost all 
galaxies possess a component of metal-poor GCs, with the average fraction of metal-rich GCs ranging 
from 15 to 60%. The colors of both subpopulations correlate with host galaxy luminosity and color, 
with the red GCs having a steeper slope. The steeper correlation seen in the mean color of the entire 
GC system is driven by the increasing fraction of metal-rich GCs for more luminous galaxies. 

To convert color to metallicity, we also introduce a preliminary (<?-z)-[Fe/H] relation calibrated to 
Galactic, M49 and M87 GCs. This relation is nonlinear with a steeper slope for [Fe/H] < —0.8. 
As a result, the mctallicities of the metal-poor and metal-rich GCs vary similarly with respect 



to galaxy luminosity and stellar mass, with relations of [Fe/H] 
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, respectively. Although these relations are shallower than the 



mass-metallicity relation predicted by wind models and observed for dwarf galaxies, they are very 
similar to the mass-metallicity relation for star forming galaxies in the same mass range. The offset 
between the two GC populations varies slowly (oc MP/ ) and is approximately 1 dex across three or- 
ders of magnitude in mass, suggesting a nearly universal amount of enrichment between the formation 
of the two populations of GCs. We also find that although the metal-rich GCs show a larger dispersion 
in color, it is the metal-poor GCs that have an equal or larger dispersion in metallicity. The similarity 
in the M*-[Fe/H] relations for the two populations, implies that the conditions of GC formation for 
metal-poor and metal-rich GCs could not have been too different. Like the color-magnitude rela- 
tion, these relations derived from globular clusters present stringent constraints on the formation and 
evolution of early-type galaxies. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: star clusters 
— globular clusters: general 
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1. INTRODUCTION 

Globular clusters (GCs) are found in nearly every 
nearby galaxy irrespective of its luminosity or gas con- 
tent. The ubiquity and relative simplicity of these old 
star clusters makes them a fundamental tool for un- 
derstanding the star formation, metal-enrichment, and 
merging histories of galaxies in the local universe. In 
particular, the color distributions of GC systems have 
played an important role in constraining the evolution 
of elliptical galaxies. While the color-magnitude dia- 
gram of early-type galaxies has garnered much interest 
as a strong constraint on the formation of ellipticals (e.g. 
Bower, Lucey, & Ellis 1992; Stanford, Eisenhardt, & 
Dickinson 1998), the mean color of an entire galaxy is 
a crude tool that necessarily combines its detailed star 
formation and chemical enrichment history into a sin- 
gle number. By comparison, globular clusters trace each 
major epoch of star formation, and because they are in- 
dividually resolved single-age, single-metallicity systems 
with typical ages older than 10 Gyr, they provide a means 
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to determine the distribution in metallicity of the early 
major star forming events that built their host galaxies. 

Metallicity distributions of the field star populations 
in most nearby ellipticals are difficult or impossible to 
obtain with current technology, therefore the color dis- 
tributions of their GC systems (inferred to be distribu- 
tions in metallicity) provide the only direct probe of their 
chemical enrichment history. With the dichotomous na- 
ture of the Milky Way's own GC system as a starting 
point (Searle & Zinn 1978), the availability of the Hub- 
ble Space Telescope (HST) over the past 15 years has 
enabled studies of the color distributions of extragalactic 
GC systems, and has revealed that bimodality is a com- 
mon property in massive elliptical galaxies (Gebhardt & 
Kissler-Patig 1999; Kundu & Whitmore 2001). This has 
led to the nomenclature of blue or "metal-poor" , and red 
or "metal-rich" GC populations, which in the Milky Way 
would correspond to the halo GCs (([Fe/H]) ~ —1.59) 
and the bulge/thick disk GCs (([Fe/H]) ~ -0.55) (Cote 
1999). Those dwarf elliptical galaxies that have been 
studied to date exhibit purely unimodal populations of 
metal-poor GCs (Lotz, Miller, & Ferguson 2004). 

Why these color distributions should be broad 
or bimodal — i.e. why there should be distinct GC 
subpopulations — is a key issue in the quest to assemble 
a consistent picture of bulge, halo, and elliptical galaxy 
formation (see West et al. 2004 for a review). At the 
very least, many if not all spheroidal systems cannot have 
formed in a single, isolated, monolithic starburst. We see 
in the local universe that massive star clusters are formed 
wherever there is a high surface density of star formation 
(Larsen & Richtler 2000), and especially in major merg- 
ers (e.g. Zhang, Fall & Whitmore 2001). Thus, some 
proposed explanations for GC color bimodality invoke 
sequential episodes of star formation to create the two 
populations, both induced by major mergers (Ashman 
& Zepf 1992) and in isolation (Forbes et al. 1997). Oth- 
ers attempt to explain GC systems in the context of 
hierarchical merging with gas dissipation in both semi- 
analytic (Beasley et al. 2002) and hydrodynamic models 
(Kravstov & Gnedin 2005). While these approaches are 
in some ways the most promising and represent the more 
quantitative work applied to this problem, neither model 
naturally produces bimodality. Cote et al. (1998, 2000) 
showed that multiple star forming events were not nec- 
essary to produce bimodal metallicity distributions in a 
hierarchical merging framework without gas dissipation. 
We see evidence of this kind of gas-poor merging both 
locally in the form of accreted dwarf galaxies like the 
Sagittarius dwarf spheroidal galaxy (Ibata, Gilmore & 
Irwin 1995), and also in mergers of massive ellipticals in 
high redshift clusters (van Dokkum et al. 1999). 

It is probable that aspects of all of these scenarios 
(which themselves are not mutually exclusive) are impor- 
tant for the formation of globular cluster systems. De- 
termining the dominant mechanisms for the formation 
of the GC subpopulations will come from quantifying 
the detailed nature of these GCs. Studies of metal-poor 
and metal-rich GCs in individual galaxies have already 
shown that they have different spatial properties with the 
metal-rich GCs being more concentrated toward the cen- 
ters of galaxies (e.g. Geisler, Lee & Kim 1996), and that 
they have different kinematic properties (Sharpies et al. 
1998; Zepf et al. 2000; Cote et al. 2001, 2003; Peng et al. 



2004). Another approach to understanding GC systems, 
which we take in this paper, is by precisely quantifying 
their variation as a function of host galaxy properties. 
In addition, this allows us to address the similarities and 
differences between GCs and field stars in a systematic 
way, providing insight into the nature of star formation 
and gas flows in the early universe. 

Previous studies have investigated GC color or metal- 
licity distributions as a function of host galaxy proper- 
ties, van den Bergh (1975) and Brodie & Huchra (1991) 
discussed the correlation between the mean metallicity of 
GC systems with the mean metallicity of their host galax- 
ies. It is now well-established that more luminous and 
metal-rich galaxies have more metal-rich GCs. Moreover, 
the mean color of the metal-rich GC subpopulation also 
appears to correlate with galaxy luminosity. Both Forbes 
et al. (1997) and Kundu & Whitmore (2001) found that 
the colors of the metal-rich GCs correlated with host 
galaxy luminosity, but neither were able to detect a cor- 
relation for the metal-poor GCs. Forbes & Forte (2001) 
found the same to be true when the GC colors were com- 
pared to galaxy velocity dispersion. 

Larsen et al. (2001) conducted a careful and homo- 
geneous study of both GC subpopulations with deep 
HST/WFPC2 photometry in the F555W and F814W 
filters of 17 galaxies, all of which were brighter than 
Mb = —18.6. In the GC color distributions of these 
galaxies, the mean colors of both the metal-rich and 
metal-poor peaks increased with host galaxy luminosity. 
They also for the first time found that the GC system 
colors had a weak correlation with host galaxy color (ex- 
pected because of the correlation with luminosity and 
the color-magnitude relation of early-type galaxies). Of 
particular interest lately is the degree to which the prop- 
erties of metal-poor GCs are correlated with their host 
galaxies. The lack or weakness of any observed trends 
have indicated to previous authors that the metal-poor 
GCs must be "universal" and have formed in similar 
sized gas fragments under similar conditions. Burgarella 
et al. (2001) compiled color data on the blue GCs for 
47 galaxies from the literature and found the trend to 
be weak or insignificant. At the other end of the galaxy 
luminosity function, Lotz et al. (2004) measured the GC 
color distributions of dwarf elliptical galaxies (dEs) using 
HST/WFPC2 imaging of 69 dEs in the Virgo and For- 
nax clusters and they too found a shallow relationship 
between GC color and galaxy luminosity. All of these 
studies relied on the best data at the time, HST/WFPC2 
imaging of galaxy samples with varying degrees of het- 
erogeneity, and almost all were based on V—I colors. The 
installation of the Advanced Camera for Surveys (ACS, 
Ford et al. 1998) now gives us the opportunity to conduct 
this sort of study in a much more precise and compre- 
hensive fashion. 

The ACS Virgo Cluster Survey (ACSVCS, Cote et al. 
2004, Paper I) is an HST/ACS imaging program of 100 
early-type galaxies in the Virgo cluster, and is designed 
for studying GC systems in a deep and homogeneous 
manner. The breadth and depth of this survey makes it 
the most complete and homogeneous study of extragalac- 
tic globular cluster populations ever undertaken. Char- 
acterizing the metallicity distributions of GC systems is 
one of the key goals of this project. Each galaxy is im- 
aged in the F475W and F850LP filters (5475 and zsso), 
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providing twice the wavelength baseline and metallicity 
sensitivity of the standard V—I color. The high spatial 
resolution and spatial sampling of the ACS Wide Field 
Camera (ACS/WFC) allows us to resolve the half-light 
radii of GCs at the distance of Virgo, facilitating the se- 
lection of GC candidates and enabling studies of the GC 
sizes themselves (Jordan et al. 2005a, Paper X). Our im- 
ages also provide high precision photometry and surface 
brightness profiles of the galaxies themselves (Ferrarese 
et al. 2005, Paper VI), and their nuclei (Cote et al. 
2005, Paper VIII), as well as distances to each galaxy 
using the method of surface brightness fluctuations (Mei 
et al. 2005a, b Papers IV and V). These data have also 
begun to show the connection between globular clusters 
and ultracompact dwarf galaxies (Ha§egan et al. 2005, 
Paper VII). Future studies (Mieske et al., in prep) will 
also use the ACSVCS data to address the topic of color- 
magnitude relations for the GCs themselves (e.g. Harris 
et al. 2005). Taken together, we are able to present a 
complete picture of these Virgo galaxies and their GC 
systems. 

2. DATA AND CATALOGS 

2.1. Data Reduction 

We have reduced our ACS/WFC images using a dedi- 
cated pipeline that is described by Jordan et al. (2004a, b, 
Papers II, III). Briefly, each image is combined and 
cleaned of cosmic rays using the Pyraf task multidrizzle 
(Koekemoer et al. 2002). We then create models of the 
galaxy in each filter for the purposes of subtracting the 
galaxy light. After model subtraction, we iterate with 
SExtractor (Bertin 1996) to mask objects, subtract any 
residual background, and do our final object detection 
using estimates of both the image noise and the noise 
present due to surface brightness fluctuations. Objects 
are only accepted if they are detected in both filters. 
After a selection on magnitude and ellipticity to reject 
obvious background galaxies, we use the program KING- 
PHOT (Jordan et al. 2005a) to measure magnitudes and 
sizes for all candidate GCs. KINGPHOT finds the best- 
fit King model parameters for each object given the point 
spread function (PSF) in that filter and at that location. 
Aperture magnitudes are measured on the best-fit PSF- 
convolved King model. For total magnitudes, we inte- 
grate the flux of the model to the limit of the PSF and 
apply a GC size-dependent aperture correction. This cor- 
rection is determined by convolving King models of dif- 
ferent half-light radii with both the Sirianni et al. PSF 
(F475W), and a PSF derived from bright stars in the 
Galactic globular cluster 47 Tuc (F850LP). For the pur- 
pose of measuring GC colors, we measure magnitudes 
within a 4-pixel radius aperture and then apply the av- 
erage aperture correction for a GC with a half-light ra- 
dius of 3 pc and a concentration of 1.5. The result is a 
catalog of total magnitudes, {g-z) colors, half-light radii 
(jh) and concentrations for each object. Magnitudes and 
colors are corrected for foreground extinction using the 
reddening maps of Schlegel, Finkbeiner & Davis (1998), 
and extinction ratios for the spectral energy distribution 
of a G2 star (Paper II). In this paper, we use g-z to mean 
5475 _ ^850 unless explicitly stated otherwise. 
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2.2. Control Fields 

One of the main problems with GC studies, particu- 
larly for fainter galaxies, is the contamination from back- 
ground galaxies. For dwarf galaxies where there may 
only be a few GCs, the contamination from the back- 
ground is often larger than the signal. For brighter galax- 
ies, the light from the galaxy itself creates a spatially 
varying detection efficiency that complicates the com- 
parison to blank fields. This is a problem that requires 
careful treatment in order to minimize systematic uncer- 
tainties, especially in low-luminosity systems. 

In order to address this issue, we use archival 
ACS/WFC imaging of 17 blank, high-latitude fields that 
have been observed with the F475W and F850LP fil- 
ters to the same or greater depth than our ACSVCS 
observations. These images were taken as part of two 
ACS Pure Parallel programs (GO-9488,GO-9575), and 
are well-suited to the purpose of understanding the back- 
ground and foreground contamination in our GC samples 
(Tabled . Each one of these images was reduced with our 
pipeline in exactly the same way as our program data. 
Since there were no large foreground galaxies in any of 
the control fields, and the exposure times were longer 
than for our ACSVCS images, all the control data go 
deeper than our program data. Hence for each galaxy, 
we created a "custom" control sample by redoing the de- 
tection as if there was a Virgo galaxy in front of of each 
blank field, and each field had the same exposure time as 
our ACSVCS fields. To do this, we used the noise images 
for our program fields that are generated in the data re- 
duction. Thus, for each control field we can produce the 
catalog of objects that we would have detected had the 
exposure times matched our ACSVCS program images, 
and had there been a particular ACSVCS galaxy in the 
foreground. 

We find that using the control fields is consistent with 
and superior to using the local background. The degree 
of field-to-field cosmic variance for GC candidates in our 
control fields is on the order of the Poisson noise. Given 
our efficient selection of GC-like objects, there are not 
many background objects that make the cut for GC se- 
lection (see below). For example, in our faintest galaxy, 
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VCC 1661, we only expect on average about ~ 9 ± 4 
contaminants to the sample of GC candidates. Deter- 
mining this background locally as opposed to over many 
control fields, while possibly more representative of the 
true background, introduces a large error due to count- 
ing noise of small numbers. In general, we find that the 
number of contaminants that we estimate from the con- 
trol fields is a good match to the number of contaminants 
found in the dwarf galaxy fields, where the number of 
objects is dominated by the background (see Peng et al. 
2005). 

2.3. Globular Cluster Selection 

The selection of a clean and complete catalog of GCs 
is an essential but difficult part of all GC studies. The 
ability to resolve and measure sizes for GCs in Virgo 
with HST provides much greater leverage to separate 
GCs from foreground stars and background galaxies. In 
addition, our custom control fields allow us to isolate 
the location of our contaminating population in the mul- 
tidimensional parameter space. This efficient selection 
allows us to avoid stringent cuts on galactocentric ra- 
dius that may strongly bias measured properties of the 
GC system. This approach to GC selection and back- 
ground estimation makes our study particularly unique 
and important for the GC systems of faint galaxies where 
the number of GCs is small compared to the number of 
background galaxies. Previous studies have all required 
a strict selection on galactocentric radius. 

While the details of the GC selection will be described 
in greater detail in another paper (Jordan et al. 2005b), 
here we briefly outline the technique. We choose a broad 
cut on color, selecting only objects with 0.5 < g-z < 2.0, 
which generously spans the age and metallicity ranges 
typical of old star clusters. The strength of our selec- 
tion lies in the size- magnitude diagram (Figure^), where 
background galaxies are typically fainter and more ex- 
tended than GCs. This figure shows the data and custom 
control fields for two galaxies in our sample. The custom 
control fields have been randomly sampled so that only 
1/17 of the objects are plotted, equivalent to what is ex- 
pected for a single field. It is already well known that 
GCs typically have a Gaussian-like luminosity function 
(Harris 1991), and that the Galactic GCs have a median 
size of r/, ~ 3 pc. This is evident in the left hand pan- 
els of Figure ^ as a population mostly separate from the 
locus of background contaminants. Notice also how the 
custom control field includes fewer objects and is shal- 
lower for the brighter galaxy. 

We apply a Gaussian kernel to the control data to pro- 
duce a nonparametric two dimensional density distribu- 
tion of contaminants. We then fit a parametric model 
to the GC locus. In magnitude, we assume a Gaussian 
distribution, and in size, we use a nonparametric ker- 
nel estimate plus a power law tail, with the two joined 
at r/j = 7 pc. This model is then fit to the data using 
maximum likelihood estimation — the fitting procedure 
allows us to account for variations in the GCLF turnover 
magnitude and differences in the size distributions from 
galaxy to galaxy. Finally, we use these two density dis- 
tributions to evaluate the probability that a given object 
is a GC. Typically, the lines of constant probability run 
diagonally from the faint, compact region of the diagram 
to the bright, extended region (seen in Figure How- 
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Fig. 1. — Globular cluster selection diagrams for four galaxies 
spanning the magnitude range of our survey. From top to bottom: 
VCC 1226 (M49, B = 9.31), VCC 2000 (B = 11.94), VCC 1528 
(B = 14.51), and VCC 1661 (B = 15.97). Left panels show the 
objects detected in the galaxy images. Right panels show a random 
selection of 1/17 of the objects detected in 17 blank sky control 
fields that have been customized for the depth reached in each 
galaxy. Notice how the control field for VCC 1226 is shallower and 
has fewer objects because the galaxy is brighter and fills the frame. 
Object colors are coded by the probability that they are globular 
clusters. Error bars represent the median photometric error at 
those magnitudes. 

ever, even at a constant probability cut at a value of 0.5, 
this dividing line is not exactly the same from galaxy to 
galaxy. In galaxies with more GCs, it is more likely that 
objects in the "contaminant" locus are actually GCs, and 
hence they will be assigned a higher probability. For the 
purposes of this paper, all of our GC catalogs will use 
objects with GC probabilities greater than 0.5. 

Recently, some investigations have revealed the ex- 
istence of faint, extended star clusters in NGC 1023, 
NGC 3384 (Larsen & Brodie 2000; Brodie & Larsen 
2002), and other nearby spiral (Chandar, Whitmore, 
& Lee 2004) and dwarf galaxies (Sharina, Puzia, & 
Marakov 2005). Because these clusters typically have 
similar sizes and magnitudes to background galaxies and 
their nature is still uncertain, our selection cuts do not 
include them in our GC sample and we do not include 
them in our analysis. Faint, extended star clusters in 
the ACSVCS galaxies will be treated in a separate paper 
(Peng et al. 2005). 

3. RESULTS 

3.1. Color Distributions of Globular Cluster Systems 

In Figure|2 we show the (g-z)o color distributions for 
the GC systems of 100 ACSVCS galaxies, ordered by de- 
creasing total B magnitude of the host galaxy. The his- 
tograms were created by binning the data with an "opti- 
mal" bin size, which is related to the width of the distri- 
bution and the total number of objects as 2IQRn~ 1 / 3 , 
where IQR is the interquartile range (the color range 



Color Distributions of Globular Cluster Systems 




o 

0.4 0.6 0. 



8 1.0 

(g 



1.2 1.4 1.6 1.8 0.40.60.8 1.0 1.2 1.4 1.6 1 

-z) (s-z) 



0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 



(.40.60.8 1.0 1.2 1.4 1.61.: 

(g- z >„ 



Fig. 2. — Globular cluster (g-z)o distributions for all 100 ACSVCS galaxies. We plot color histograms of GC candidates (black), expected 
contaminants (cyan), and statistically cleaned GC distributions (green). The black curve represent a nonparametric kernel density estimate 
of the cleaned distribution. In cases where the distribution is likely bimodal, we plot the red and blue Gaussian components as determined 
by the KMM estimates, as well as their sum (purple). In the unimodal cases, we plot the best fit single Gaussian to the entire distribution 
(purple). 



that includes the second and third quartile of the ranked 
color distribution), and n is the total number of objects 
(Izenman 1991). The bin size is not allowed to be smaller 
than the mean photometric error. The histograms of the 
GC data are shown in black, and histograms of the ex- 
pected contamination as measured from the custom con- 
trol fields are shown in cyan. 

While the level of contamination is negligible for the 
brighter galaxies, the background is a significant problem 
for the fainter galaxies. We create statistically cleaned 



samples of GCs by using a Monte Carlo procedure. For 
each GC, we calculate the probability that it is a con- 
taminant by using a nonparametric density estimate of 
the control data as compared to the program galaxy data 
at that color. Based on this probability, we randomly in- 
clude or do not include this object (with replacement) 
from our generated sample. Iterating 100 times for each 
galaxy, we can produce an average color distribution that 
is cleaned of contaminants. These histograms are plot- 
ted in green. Kernel density estimates of the cleaned 
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Fig. 2. — continued. Globular cluster (g-z)o distributions. 

distribution are overplotted as black curves. 

The cleaned and smoothed color distributions show 
obvious bimodality for most of the brighter galaxies, 
something which is expected from previous investigations 
(e.g. Gebhardt & Kissler-Patig 1998; Larsen et al. 2001, 
Kundu & Whitmore 2001). For intermediate luminosi- 
ties and fainter galaxies, the situation becomes less ob- 
vious as the number of GCs per galaxy decreases. Yet, 
it is clear that the fraction of red GCs decreases as the 
galaxy itself becomes fainter. 

The properties of the GC color distribution are of- 
ten quantified using the Kaye's Mixture Model (KMM; 



McLachlan & Basford 1988; Ashman, Bird, & Zepf 1994) 
to fit two Gaussians to the data using the expectation- 
maximization (EM) method. We fit the homoscedastic 
case of this model, where a is the same for both Gaus- 
sians. Constraining a reduces the number of free param- 
eters and gain leverage on noisy data in which there are 
generally large errors when fitting for a individually. We 
apply KMM to the each of the probabilistically cleaned 
GC color samples to determine the means of the blue and 
red peaks in the color distribution. We determine that 
the two Gaussian model is a better fit to the data than 
the one Gaussian model if the "p-value" for the bimodal 
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Fig. 2. — continued. Globular cluster (g-z)o distributions. 

model is less than 0.05. For the cases where bimodal- 
ity is deemed significant (p < 0.05), Figure shows the 
two Gaussians and their sum. For the unimodal case, 
the single best-fit Gaussian is plotted. In some galaxies, 
the distribution was determined to be bimodal, but the 
number of GCs in the red peak was not significant over 
the background at greater than 2a. In these cases, we 
reject the bimodal hypothesis because the red objects are 
likely to be background galaxies. The results for KMM, 
and the parameters for the best fit one or two Gaussian 
models are presented in Table El This table includes the 
means of the blue and red peaks, the common sigma, 



the fraction of GCs that are determined to be in the red 
subcomponent, the total number of GCs (after account- 
ing for the background), the p- value of the two Gaussian 
hypothesis and its associated error, and the mean and 
sigma for the entire distribution. In cases where the red 
peak was less than 2a above the background, only the 
mean color of the blue peak is reported. In cases where 
the total number of GCs was statistically equal to zero 
(see next section), no values could be reliably estimated. 



3.2. Notes on Special Galaxies 
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Fig. 2. — continued. Globular cluster (5-2)0 distributions. 

A few galaxies are noteworthy in that they appear 
somewhat different from others of similar luminosity. For 
a more in-depth description of every galaxy in the sam- 
ple, see Ferrarese et al. (2005, Paper VI). Four galaxies in 
particular, despite their low luminosities, have prominent 
blue and red peaks and have large numbers of GC. This 
is due to their proximity to the two giants of the clus- 
ter, M87 and M49, and we are likely observing the GC 
systems of their larger neighbors. Interestingly, all four 
are also significantly redder than one would expect for 
galaxies of their luminosity, and both the number and 
mean color of the red GCs in these galaxies is similar 




0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

(g-z) (g-z) () 



to those seen in galaxies of the same (g-z) color. This 
suggests that they may be remnants of what were once 
larger, more luminous systems. 

VCC1327/NGC4486A: This dwarf elliptical is only 
7.5' away from M87 (NGC 4486), and its GCs are likely 
to be dominated by those of the nearby cD galaxy. 
VCC 1297/NGC4486B: This compact elliptical 
galaxy is similar to M32 in appearance, and is only 7.3' 
away from M87. Some or all of its GCs may in fact 
belong to the M87 GC system, and its current size and 
luminosity may be the result of significant tidal stripping. 
VCC 1192/NGC4467: This elliptical galaxy is also 
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Fig. 2. — continued. Globular cluster (g-z)o distributions. 

compact in appearance and is only 4.2' away from M49 
(NGC 4472). Its GC system may be dominated by that 
of M49. 

VCC 1199: This elliptical is also only 4.5' from M49. 
Like the other small galaxies near giant ellipticals, it dis- 
plays a prominent red peak despite its low luminosity. 
Other galaxies of note are: 

VCC 1146/NGC4458: This galaxy is unique in the 
sample in that its GC system is dominated by a single 
red peak of GCs. The red GC fraction of 0.84 estimated 
by KMM is not only much larger than the value of 0.3 
expected for a galaxy of this luminosity, but is much 



higher than the 0.6 fraction seen in the giant ellipticals. 
VCC 798: This galaxy appears to be the best candidate 
in our sample for having a trimodal color distribution. 
VCC 731: Appears to have significantly more GCs than 
other galaxies of comparable luminosity. The excess ap- 
pears to be due to an large number of red GCs. 
VCC 1692: The blue and red GCs are particularly well 
separated in color. 

VCC 1499: While this galaxy appears to be an elliptical 
on ground based Digitized Sky Survey images, our ACS 
images reveal it to contain numerous young blue star 
clusters at the center. Also, the color of the galaxy is 
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TABLE 2 

Color Distributions: Bimodal and Unimodal Parameters 



No. 
(1) 


VCC 
(2) 


f^blue 

(3) 


ftblue.err 

(4) 


ftred 

(5) 


/^red.err 

(6) 


1 


1226 


0.97 


0.01 


1.42 


0.01 


2 


1316 


0.98 


0.01 


1.43 


0.01 


3 
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0.98 


0.01 


1.45 


0.01 


4 
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0.98 


0.03 


1.33 


0.03 


5 
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1.02 


0.02 


1.34 


0.03 


6 
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0.97 


0.01 


1.36 


0.01 


7 
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0.98 


0.01 


1.36 


0.01 


8 
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0.94 


0.01 


1.41 


0.01 


9 
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0.94 


0.02 


1.33 


0.01 


10 
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0.01 


1.39 


0.01 
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0.96 


0.02 


1.35 


0.01 
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0.92 


0.05 


1.22 


0.04 
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0.03 
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0.04 
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0.04 
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0.03 






24 
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0.04 
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TABLE 2 

Color Distributions: Bimodal and Unimodal Parameters 
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Fig. 3. — Mean colors of the entire GC system and the two GC subpopulations as a function of Mb- We plot the values of the blue 
and red means as determined by the KMM algorithm assuming a model of two Gaussians with a common dispersion. Only 54 galaxies 
for which the red peak was more than 2-sigma above the background are plotted with two peaks (circles and diamonds). Filled symbols 
represent the 24 galaxies that have significant bimodal distributions (p < 0.05), and the solid lines are fits to these points. Open symbols 
are for the 30 galaxies with 0.05 < p < 1 and the dotted lines are fit to data of all p-values. The means of the entire GC color distribution 
are also plotted (green and gray triangles) with 46 of the points (green) representing the GC systems determined to be unimodal in color. 



one of the bluest in the sample, suggesting that it may 
be a dl/clE transition object. In the same field, only 
1.4' away, is a true dE, VCC 1491, hence the clusters 
measured in this field are a combination of the systems 
of these two galaxies. 

VCC 9: This galaxy has a very low surface brightness 
for its reported luminosity. 

VCC 1938: This SO galaxy is 1.7' from the neighboring 
nucleated dE, VCC 1945. We detect GCs in both galax- 
ies, although there are many fewer associated with VCC 
1945 as it is ~ 2.5 mag fainter. 

VCC 33, 1779, 1627, 1993: In the full field data for 
these galaxies, the number of GC candidates we detect 
is less than 3a above the expected background. VCC 
1779 shows spiral structure and dust, is likely to have a 
younger age, and there are a two likely star clusters near 
the center. The other three galaxies do not show any 
obvious concentrations of GCs at their centers although 
there may be one or two GCs in each galaxy. Only for 
VCC 1627 do we not detect any GC candidates above 



the background. 

For the purposes of the analysis that follow, we exclude 
from our sample the galaxies VCC 1327, 1297, 1192, 1199 
(near giants), 1499, 1779 (younger ages), and 1938 (two 
galaxy blend) because either the detected globular clus- 
ter systems may not be representative of the targeted 
galaxy, or the young ages belie the classification as an 
early-type galaxy. This leaves us with 93 galaxies in our 
sample. 

3.3. GC colors and Host Galaxy Properties 

3.3.1. Color Decompositions of Individual GC Systems 

Using the mixture model estimates presented in Ta- 
ble[5J we can investigate the behavior of the blue and red 
GC subpopulations as a function of galaxy luminosity. 
Previous studies along these lines were mainly based on 
HST/WFPC2 data (Forbes et al. 1997, Kundu & Whit- 
more 2001; Larsen et al. 2001; Burgarella et al. 2001) and 
included very few galaxies with Mb fainter than —19.5. 
The study of Lotz et al. (2004) targeted dEs in Virgo 
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Fig. 4. — This image shows the GC color distributions of the ACSVCS galaxies ordered by host galaxy Mb- Each column is a kernel 
density estimation of a galaxy's GC color distribution, with white normalized to represent the peak density. The galaxies are ordered in 
increasing luminosity from left to right. The image shows that all galaxies seem to possess a population of blue GCs with similar color, 
but that the color and relative fractions of red GCs is a strong function of galaxy luminosity. 

We are able to resolve the GC subpopulations for 
galaxies ~ 2 mag fainter than those in previous stud- 
ies, and we observe a clear trend that the mean colors of 
both blue and red GCs are redder for more luminous host 
galaxies, and that the slope of the relation for red GCs 
is steeper by a factor of 1.6-1.9. We fit these relations to 
both the p < 0.05 sample and the full sample, deriving 
the following weighted linear fits of the form 



and Fornax, and saw no evidence for bimodality in their 
GC color distributions. The higher metallicity sensitiv- 
ity of the {g—z) color and the deeper photometry of our 
ACS observations makes this an ideal program for test- 
ing these relationships in a homogeneous fashion across 
a wide range of galaxy luminosities and colors. 

In Figure 01 we show the mean colors of the blue and 
red GC subpopulations as a function of the absolute 
blue magnitude of their host galaxy. The galaxy mag- 
nitudes are from the photometry of Binggeli, Sandage & 
Tammann (1985) and are listed for the entire sample in 
Paper I. We assume a distance to the Virgo cluster of 
16.5 Mpc (Tonry et al. 2001) with a distance modulus of 
31.09 ± 0.03 mag from Tonry et al. (2001), corrected by 
the final results of the Key Project distances (Freedman 
et al. 2001; see also discussion in Mei et al. 2005b). The 
GC color data follow Table and the two individual 
means are plotted for all galaxies that had a significant 
(> 2a) red GC population. For 24 galaxies plotted as 
solid points, the distributions were significantly bimodal 
(p < 0.05). The color distributions for 30 other galaxies 
(open points) can also be decomposed into two compo- 
nents, but are less uniquely described by a two-Gaussian 
model (0.05 < p < 1). In both cases, we plot the means 
of the two fitted Gaussians only if the number of red 
GCs is more then 2a above the expected background. 
This eliminates small numbers of red background galax- 
ies from causing spurious bimodality. For the remaining 
46 galaxies (those that have insignificant numbers of red 
GCs), we treat them as unimodal and plot the means of 
their entire GC color distributions (triangles). 



(g — z) = a + b x Me 



(1) 



Coefficients for these fits are listed in Table [3J While 
most galaxies with Mb < — 18 have distinguishable blue 
and red components in their GC color distributions, the 
fraction of galaxies whose distributions are not well re- 
solved increases for fainter galaxies until all of the galax- 
ies with Mb > — 16 have an insignificant number of red 
GCs. Because both the fraction of red GCs and the to- 
tal number of GCs is lower for fainter galaxies, the red 
subpopulation becomes more difficult to separate out at 
intermediate luminosities. However, the slopes that we 
derive from the brighter galaxies are statistically iden- 
tical to those fit to the larger sample. In other words, 
even when a galaxy's GC color distribution is not sta- 
tistically very different from unimodal, a decomposition 
into two components is still consistent with the trends 
seen in more luminous galaxies, and thus may still be an 
appropriate description of the system. 

3.3.2. Color Distributions Binned by Galaxy Properties 

When the data is of high enough quality, it is some- 
times best to view the data itself in aggregate rather than 
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Fig. 5. — Kernel density color distributions of GCs in seven bins of host galaxy magnitude (a) and six bins of host galaxy (g— z)q (b). 
Magnitude bins are 1 mag wide and extend from Mg = —21 to Mb = —15 (red to blue). The distance to the Virgo cluster is assumed to 
be 16.5 Mpc. Bins of color are 0.1 mag wide from 1.6 to 1.1, with the last bin extending to 0.75 (red to blue). Distributions have been 
corrected for the expected background and foreground contamination as measured from our control fields. All color distributions have 
been constructed with a Gaussian kernel (a = 0.05 mag). Notice how even the faintest and bluest bins in our sample appear asymmetric, 
indicating the presence of some red GCs. Data for this plot is listed in Tables HI and 151 



struggle with parametrization. Figure 0] displays a non- 
parametric representation of the GC color distributions 
for 92 galaxies (the 93 from our analysis sample minus 
VCC 1627, which has no detected GCs). Each column of 
this plot is the background cleaned color distribution of a 
single galaxy's GCs constructed with a Gaussian kernel, 
where the grayscale has been scaled such that zero den- 
sity is black and the mode of the distribution is white. 
The galaxies are rank ordered by their luminosity with 
approximate Mg labeled on the x-axis. Immediately ap- 
parent from this image are the two GC subpopulations 
and their behavior with galaxy luminosity. Nearly all 
galaxies appear to possess a blue subpopulation of GCs, 
and the mean of this population varies only slowly with 
galaxy luminosity. In addition, there is a population of 
red GCs whose color and number fraction increase in 
a continuous fashion across our entire sample, spanning 
seven magnitudes of galaxy luminosity. Lotz et al. (2004) 
made a similar inference when comparing the GCs in dEs 
to the blue GCs in ellipticals. Interestingly, a red wing 
in the GC color distribution appears to exist even in 
some of our faintest galaxies, although the small number 
of GCs and increased noise due to background subtrac- 
tion makes it difficult to quantify this for any individual 
galaxy. 

Figure 0] displays the utility of grouping galaxies to- 
gether by their intrinsic properties. Because of the large 
number of galaxies in our sample, we can quantify trends 
by treating the galaxies collectively even when any sin- 
gle faint galaxy has too few GCs for substantive analysis. 
In Figure [3 we investigate the GC color distributions as 
a function of both galaxy luminosity and galaxy color 
by binning them and accumulating enough GCs to over- 
come the noise. We create seven bins of magnitude from 
M B = —22 to —15 in steps of one magnitude. We have 
six bins in galaxy color, five from (g — z)q = 1.6 to 1.1 



in steps of 0.1 mag, and the last bin somewhat larger 
0.75 < (g — z)o < 1.1. The colors of the galaxies were 
derived from our ACS images as described in Paper VI. 
Briefly, the colors were estimated by directly integrating 
the observed surface brightness profiles between 1" and 
the smaller of one effective radius or the radius at which 
the surface brightness profile falls one magnitude below 
the sky (in either filter). 

The cumulative background cleaned GC color distri- 
butions for these galaxy bins are presented in Figure [SJ 
The trends seen in Figure 0] are also evident in the first 
plot of Figure El In each color distribution a strong blue 
peak is easily visible, and each distribution is asymmetric 
about this blue peak, even the one for our faintest bin. 
The behavior of the red peak is also easy to discern. The 
number fraction and color of the red peak progressively 
decreases for fainter galaxies. Where the red peak is ob- 
vious, it also appears to be broader than the blue peak. 
The same progression is seen for bins in galaxy color. 
Redder galaxies have more red GCs, and these red GCs 
are themselves redder. 

Each of the magnitude and color bins across our entire 
sample possess a population of blue GCs with some addi- 
tional red GCs. We can try to decompose these two pop- 
ulations in our binned sample in a way similar to how we 
treated the individual galaxies, by treating them as the 
sum of two Gaussian distributions. However, given the 
higher signal-to-noise of these summed distribution, we 
find that the sum of two Gaussians, whether with identi- 
cal or independent variances, is often not a good model 
for the color distributions for the purpose of measuring 
the colors of the two modes. When we run the KMM al- 
gorithm on these data to test for bimodality, all the color 
distributions return small p- values, including the faintest 
and bluest bins. However, when estimating the mean 
values of the blue and red GC peaks, the two-Gaussian 
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Fig. 6. — Color histograms and nonparametric decompositions 
of GCs binned by host galaxy magnitude. 

fits systematically skew the colors of the blue and red 
peaks toward the median of the entire distribution — i.e. 
the estimated color of the blue peak is too red, and the 
estimated color of the red peak is too blue as compared 
to a non-parametric estimation from the data. This ef- 
fect is most pronounced when the modes are of nearly 
equal strength, and the result is to introduce a bias that 
is correlated with galaxy luminosity or color. This is 
something that needs to be treated with care because it 
can produce correlations that are merely artifacts of the 
parameter estimation. 

3.3.3. Nonparametric Decompositions of Binned GC 
Color Distributions 

While the addition of another Gaussian component 
(and its attending free parameters) does improve the fit 
with marginally significant p- values in some cases, we 
instead choose to avoid this somewhat arbitrary model 
and decompose the two peaks in a nonparametric way. 
Given the ubiquity of the blue GCs, we start with the 
assumption that every color distribution is made up of 
a population of blue GCs with a symmetric distribution 
about the blue mode. We then take the GCs blueward 
of the blue mode, reflect them about the blue mode, and 
take that to be the blue GC population. The remain- 
der of the GCs are then what we consider the red GCs. 
The peak color of the red peak is estimated from the 
GC color distribution that has the blue GCs subtracted. 
While this method does make assumptions, particularly 
about the symmetry of the blue GC color distribution 
and the choice of a kernel size, unlike a multi-Gaussian 
model it makes no assumptions about the shape of the 
red GC color distribution. The results are also not par- 
ticularly sensitive to the size of the kernel as long as it is 
not so small as to introduce spurious peaks, and not so 
large as to be comparable to the half-width of the blue 
peak. 

In Figures IHI and we overplot our nonparametric de- 
compositions for the GC color distributions, binned by 
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Fig. 7. — Color histograms and nonparametric decompositions 
of GCs binned by host galaxy color. 



galaxy magnitude and color. We estimate the median 
color and the half-width that encompasses 68% of the 
GCs for each peak, and also for the entire GC popula- 
tion in that bin. The parameters and their errors for 
this model were estimated using the bootstrap with 1000 
iterations. These values are listed in Tables 0] and [S] 

Using these decompositions, we can study the proper- 
ties of the GC color distribution as a function of galaxy 
properties. The results of these fits are shown in Fig- 
ures !8ll0l In Figure|S|i we show the relationship between 
the median colors of the red, blue, and total GC popula- 
tions against galaxy luminosity and galaxy color. As in 
FigureOl the colors of the red, blue, and total GC popu- 
lations correlate with the luminosity of the host galaxy. 
The coefficients for the linear fits are presented in Ta- 
ble El The relationship versus luminosity is similar to 
that shown with individual galaxies in Figure [3J How- 
ever, there are two notable differences. First, we are able 
to trace the mean colors of the red and blue peaks for 
the entire magnitude range of our sample with much less 
noise. Second, the slopes of the red and blue relations 
are more disparate than they are in Table The slope 
for the red GCs is steeper and the slope for the blue GCs 
is shallower than was derived from individual galaxies. 
The result is that the slope for the red GCs is 4.6 times 
steeper than for the blue GCs. 

Why would the two different methods give different 
values? This difference stems from two effects. The first 
is that the individual galaxy and binned galaxy data were 
decomposed using two different methods. When we ap- 
ply the KMM algorithm to the binned data, we do get 
a steeper blue GC relation and a shallower red GC rela- 
tion. However, as we noted before, this is at least in part 
due to biases in the fitting of the peaks that arises from 
non-Gaussianity in the underlying distributions. While 
a homoscedastic two-Gaussian fit is often the best one 
can do for the GCS of a single galaxy, the binned data 
provide a more critical test of an inadequate model, and 
also make nonparametric methods more feasible. 
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The second reason for the difference is that the GC 
systems of individual galaxies sometimes cannot be de- 
composed into two populations, either because they have 
no red GCs, the number of red GCs has low signifi- 
cance, or the total number of GCs is simply insufficient. 
These galaxies are necessarily dropped, thus biasing the 
data toward galaxies with well-spaced GC subpopula- 
tions. This is a bias that will be strongest at the faint 
end of the galaxy luminosity function, and will result in 
a shallower slope for the red GCs and a steeper slope 
for the blue GCs. The binned data includes all galaxies, 
regardless of their individual decomposition, and thus is 
less biased in this fashion. 

Figure QJ) also shows the relationship between the col- 
ors of these peaks with the color of the galaxy. The bluest 
bin is somewhat suspect as it contains few galaxies and 
is quite noisy. The other bins, however, each contain 
> 400 GCs. Both of the GC subpopulations show a clear 
correlation with galaxy color, with the red GCs showing 
a steeper relation. Larsen et al. (2001) found that the 
two peaks correlate with galaxy color at the 2-4<7 level 
(see their note added in manuscript), and the higher pre- 
cision of the ACS photometry for both the GCs and the 
galaxies leaves little to doubt. While this relationship 
might be expected since galaxy color is known to cor- 
relate with luminosity, the color-magnitude relationship 
for galaxies has significant scatter at the faint end (see 
Ferrarese et al. 2005) and it is not obvious that a tight 
correlation should necessarily exist. The coefficients of 
the weighted linear fits are given in Table The dashed 
line in this plot shows the one-to-one relation between 
host galaxy and GC colors. For bright galaxies, the me- 
dian color of the entire GC system appears to track the 
galaxy color closely with a ~ 0.3 mag offset, but the slope 
changes for faint galaxies. 

In Figure El we show the fitted color dispersion in each 
subpopulation, and the dispersion of the entire color dis- 
tribution as a function of galaxy luminosity and galaxy 
color. In both plots, we can see that the dispersion 
in the GC colors increases for more luminous and red- 
der galaxies. Although most work on individual galaxies 
necessarily assume that the widths of the blue and red 
distributions are the same, these decompositions for our 
binned distributions shows that, at least for the brighter 
galaxies, the dispersion in color of the two populations 
may not in fact be the same but that the red GCs may 
have a larger spread in color. We note that the way we 
have chosen to do our nonparametric decomposition, the 
width of the blue GC distribution is entirely dependent 
on the half of that is blueward of the peak. However, 
KMM estimates of the dispersions using a heteroscedas- 
tic Gaussian model give a similar result. 

One of the other noticeable trends across galaxy lumi- 
nosity and color involves the fraction of red GCs. Fig- 
ure shows how the fraction of red GCs is much higher 
in more luminous and redder galaxies. We quantify this 
in Figure 1101 showing that even in the faintest bin in 
our sample, galaxies on average have a ~ 15% fraction 
of red GCs, and this increases to ~ 60% for the most 
luminous galaxies. In the most luminous galaxies, how- 
ever, we are only sampling the inner regions of the GC 
system and color gradients could slightly affect the total 
fractions (see e.g. Rhode & Zepf (2004) for a wide-field 
study). Nevertheless, This changing fraction of red GCs 
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Fig. 8. — Mean values of blue, red, and total GC colors as a 
function of host galaxy luminosity (left) and color (right). Magni- 
tude and color bins, are the same as in Figure [F] Points represent 
the blue peak of the GC color distribution (blue diamonds), the 
red component (red filled circles), and the mean of the entire dis- 
tribution (asterisks). The dotted lines are linear fits to the two 
subpopulations, and the dashed line (right) has slope unity and 
represents the galaxy colors. Both populations of GCs get redder 
with galaxy luminosity and color, but the slope for the red GCs is 
many times steeper. 
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Fig. 9. — Color dispersions (s68, the 68% half-width) of blue 
(diamond), red (dot), and total (asterisk) GC distributions as a 
function of host galaxy luminosity (left) and color (right) binned 
by magnitude. Magnitude, and color bins same as in Figure [5] 
For luminous and red galaxies, the width of the red subpopulation 
appears to be larger that that for the blue GCs. 



Color Distributions of Globular Cluster Systems 



17 



1.0 



0.8 - 



.2 0.6 - 



U 

a 



0.4 - 



0.2 



0.0 



-16 -17 



-18 -19 
M 



-20 -21 



1.00 1.10 1.20 1.30 1.40 1.50 
Host Galaxy (g-z) 



Fig. 10. — Fraction of red GCs as a function of host galaxy lu- 
minosity, binned by magnitude. The fraction of red GCs increases 
by over a factor of six across the magnitude range of our sample. 

is the main driver behind the correlation between the 
mean color for the entire system and the luminosity of 
host galaxy. 

4. DISCUSSION 
4.1. The (g-z)-[Fe/H] Relation for Globular Clusters 

In order to make inferences on the nature of GCs and 
their host galaxies, it is necessary to translate observables 
(color, magnitude) into physical properties (metallicity, 
luminosity, mass). One of the more critical transforma- 
tions is that of color to metallicity for globular clusters. 
Because GCs are old, simple stellar populations, broad- 
band color should be a good proxy for metallicity. Be- 
yond an age of a few Gyr, the color of a star cluster 
is predominantly determined by its metal content, and 
because it was formed in a single burst, there is no com- 
plicated star formation history to disentangle. However, 
in practice, determining the precise transformation from 
different filter systems to [Fe/H] even for this simplest of 
stellar populations is a difficult task. 

Evolutionary synthesis models of stellar populations 
can in theory produce these relations, but broadband 
colors provide a challenge as they are highly dependent 
on the spectrophotometry of the input stellar libraries, 
whether theoretical or empirical. Recent models (e.g. 
Bruzual & Chariot 2003) show reasonable agreement 
with Milky Way, M31, and Magellanic Cloud star clus- 
ters. Others have taken an empirical approach and fit 
various relationships to the measured colors and metallic- 
ities of GCs in the Milky Way and other nearby systems. 
The color used most often is V— I and numerous relations 
have been calibrated, although their slopes can vary by a 
factor of two (Couture, Harris, & Allwright 1990; Kisslcr- 
Patig et al. 1997; Kissler-Patig et al. 1998; Kundu & 
Whitmore 1998; Barmby et al. 2000). All of these cali- 
brations rely heavily on the Milky Way GCs (although 
Kissler-Patig 1998 includes GCs from NGC 1399 to sup- 
plement the metal-rich end). The main limitations of this 
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Fig. 11. — [Fc/H] versus (g-z)o for 40 low-extinction Milky Way 
GCs, 33 M49 and M87 GCs with ACSVCS photometry, and 22 
M49 and M87 GCs with SDSS photometry. All GCs have published 
spectroscopic metallicity determinations as described in the text. 
The relatively low scatter in the data compared to previous work 
shows that the relation is very likely non-linear, with a steepening 
at [Fe/H] < —0.8. The dotted line is a broken linear bisector fit to 
all GCs with [Fe/H] < 0. 

approach are the homogeneity of the Galactic GC pho- 
tometry the accuracy of the published reddenings, and 
the inherent sensitivity of the V-I color to metallicity. 
Moreover, models predict that the color-metallicity re- 
lationship should be nonlinear, especially for metal-poor 
populations where metallicity sensitivity in the optical is 
diminished. 

The color that we are working with, (g-z), is twice as 
sensitive to metallicity as V-I (Paper I), and thus offers 
higher metallicity resolution than previous studies. Be- 
cause of the importance of deriving this transformation, 
we have completed a program to image the Milky Way 
GC system in these two bandpasses in order to derive an 
empirical (g-2;)-[Fe/H] relation. We used the Cassegrain 
Focus CCD Imager on the CTIO 0.9-meter telescope over 
two observing runs (8 May - 12 May 2003 and 31 May - 
6 June 2004) to observe ~ 100 Milky Way GCs in the g' 
and z' filters. We present a preliminary color-metallicity 
relation in this section, while the details and final analy- 
sis of this data set will be presented in a separate paper 
(West et al., in prep). 

From a sample of 76 Milky Way GCs with good pho- 
tometry, we selected 40 that have low reddenings of 
E(B - V) < 0.3. Reddenings and [Fe/H] values for 
these GCs were obtained from the McMaster Milky Way 
GC catalog (Harris 1996) with reddenings and metallici- 
ties predominantly compiled from Reed, Hesser, & Shawl 
(1985), Webbink (1985), Armandroff & Zinn (1988), and 
Zinn (1985). To supplement this sample, especially at 
higher metallicities, we added GCs in the giant ellipti- 
cals M87 and M49 that have spectroscopic metallicities 
(Cohen, Blakeslee, & Ryzhov 1998; Cohen, Blakeslee, & 
Cote 2003). The latter metallicities were measured using 
models of Worthey (1994) and calibrated to the metal- 
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licity scale of Zinn & West (1984). For 33 of these GCs, 
we have (3475-^850) colors from our ACS/WFC photom- 
etry. For another 22, we were able to retrieve gsDSS and 
zsdss photometry from the Sloan Digital Sky Survey 
Third Data Release (Abazajian et al. 2005), and cor- 
rect them for reddening from the maps of Schlegel et al. 
(1998). All photometry was shifted to the HST/WFC 
AB photometric system (3475 and zsso)- 

In Figure El we show the combined sample of 95 glob- 
ular clusters with both a metallicity determination and 
a (g-z) color. While there is still sizable scatter in the 
data, the relation is tight enough to show a significant 
departure from linearity with the relation steepening for 
[Fe/H]< —1.0. The dotted line shows a broken linear fit 
to the data, with a break point of (g— z)= 1.05. Because 
there is appreciable scatter in both axes, we fit the rela- 
tion as both [Fe/H] vs. (g-z) and (g-z) vs. [Fe/H], and 
take the mean of the results (e.g. Barmby et al. 2000). 
We take this approach because there is both intrinsic and 
observational scatter in both variables, and the scatter 
is not well-quantified. We also fit the data using an or- 
dinary least squares bisector fit and the method of bi- 
variate correlated errors with intrinsic scatter (BCES; 
Akritas & Bershady 1996). We find that both of these 
methods give similar results, with the slopes varying by 
less than 4% depending on the chosen method, a range 
much smaller than the error (see also Isobe et al. 1990 
for a discussion of linear regression methods). Our fit 
does not include the GCs with [Fe/H]> 0. Not only are 
the Worthey (1994) models not well-calibrated in that 
regime, but the mean colors of the red GCs that we are 
concerned with in this paper always have (g—z)o < 1.5. 

The relation we derive is: 

IFp /"Hi — / ~ 6 - 21 + (5 14 ± °- 67 )<S475 - 2=850). 0.70 < (9475 - *85o) < 1-05 
[re/ XXj — ^ _ 2 75 + (1 83 ± . 23)(s475 _ , 85o)i 1.05 < (g475 _ Zfj50 ) < 1.45 

(2) 

In Figure 1121 we show the same data and the broken 
linear fit, but plot predicted colors from three different 
models with different initial mass functions (IMF). We 
choose the high resolution models of Bruzual & Char- 
lot (2003) with IMFs of Chabrier (2003) and Salpeter 
(1955), the 2001 release of the Bruzual & Chariot (1993) 
models, and the PEGASE models, Version 2 (Fioc & 
Rocca-Volmerange 1997) with both Salpeter and Kroupa 
(2001) IMFs. In all cases, we have linearly interpolated 
the model data (black points) for a 13 Gyr old simple 
stellar population. While there is decent agreement be- 
tween the models at metallicities [Fe/H]> —1, the rela- 
tions are increasingly in disagreement at lower metallici- 
ties, and the steepness of the relation makes the choice of 
IMF an important one. Additionally, the lack of model 
isochrones for —1.6 < [Fe/H] < —0.7 means that there 
is little information on the color-metallicity relation for 
a range that is crucial for GCs, and where the slope of 
the relation is rapidly changing. Because of this, we will 
adopt our preliminary empirical relation to obtain metal- 
licities in this paper. Nevertheless, we urge caution to all 
who would transform GC colors into metallicities using 
either theoretical or empirical relations. 

Lastly, we note that by using the Milky Way GCs to 
calibrate the color-metallicity relation, we are assuming 
that all globular clusters have old (> 10 Gyr) ages. This, 
however, may not be the case in all GC systems. For 
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Fig. 12.— [Fe/Hl versus (g-z)o for Milky Way, M49, and M87 
GCs (as in Figure [TO with stellar population models. For each 
model we plot the [Fe/H] and (g-z)o values for a 13 Gyr simple 
stellar population (black points, linearly interpolated) . The solid 
purple line is the broken linear bisector fit from Figure ITT1 While 
the models show decent agreement with the data, there is a large 
spread at low metallicity, and there are also no models for metal- 
licities in the crucial range —0.6 < [Fe/H] < —1.6. 

example, spectroscopic age determinations of metal-rich 
GCs in a few other galaxies show that some of them 
can have intermediate ages of 3-8 Gyr (Goudfrooij et al. 
2001; Peng et al. 2004; Puzia et al. 2005). For the trends 
we see in color to be caused solely by age trends, though, 
would require that the red GCs in our faintest galaxies 
have a mean age of 3 Gyr, assuming that the GCs in the 
most massive galaxies were 13 Gyr old. While this mean 
age is unlikely given the current data, it is not out of 
the question that smaller galaxies may have somewhat 
younger metal-rich GCs. The age difference necessary 
for the blue GCs would be smaller (7 Gyr old in our 
faintest bin). However, all spectroscopic evidence to date 
indicates that metal-poor GCs in nearby galaxies have 
ages that are consistent with those of the Galactic GC 
system. Thus, we will assume that the Milky Way GCs 
can provide a good calibration for extragalactic GCs with 
the caveat that if there is a trend for less luminous or 
bluer galaxies host younger GCs, it will serve to flatten 
the relations we derive for metallicity. 

4.2. Globular Cluster Metallicities, Galaxy Luminosity, 
and Stellar Mass 

Using our empirical color-metallicity relation, we can 
transform the colors of the binned GC populations to 
[Fe/H]. In Figure llTfl we plot the metallicities of the blue, 
red, and total GC populations against the host galaxy 
Mb- The coefficients for the linear fits to these relations 
are presented in TableEJ As was suggested by Figure|BH, 
the metallicities of both the metal-rich and metal-poor 
GC populations increase for more luminous host galax- 
ies. One noticeable difference from Figure |Hk, however, 
is that the slope of the relation for the metal-poor GCs 
is more pronounced, and is similar to the slope for the 
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Fig. 13. — Mean metallicities of blue, red, and total GC pop- 
ulations as a function of galaxy luminosity. Bins and points are 
the same as for Figure [S] The asterisk represents the spheroidal 
component of the Milky Way (bulge and halo) and its GC system. 
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Fig. 14. — Mean metallicities of blue, red, and total GC popu- 
lations as a function of galaxy stellar mass. Stellar mass has been 
translated from luminosity using an average, luminosity-weighted 
mass-to-light ratio. Data for this plot is in Table 171 



metal-rich GCs. Given the (g-z)-[Fe/H] relation in equa- 
tion |21 the mctallicity of the GG populations are pro- 
portional to a power of the luminosity as Z oc L a where 
a = 0.16±0.04, 0.26±0.03, and 0.52±0.02 for the metal- 
poor, metal-rich, and total GC populations. These errors 
include the formal errors in the color-metallicity relation. 
Lotz et al. (2004) fit the relation to the V-I colors of 
the total GC populations in Virgo and Fornax dEs (as- 
suming only a single component), and their slopes were 
Z oc £0.i6±o.05 tQ L o.22±o.05 depending on whether they 
included Local Group dwarf spheroidals or the blue GCs 
from Larsen et al. (2001). This is in reasonable agree- 
ment with our fits especially considering that our sample 
galaxies are generally more massive, and that they used 
both a different color and color-metallicity relation. 

One of the most fundamental galaxy properties is its 
mass. While we do not have dynamical masses for these 
galaxies yet, we can obtain a rough estimate of its stellar 
mass-to-light ratio from its broadband optical-infrared 
colors. We use our (g—z) colors and J-K from the Two 
Micron All Sky Survey Extended Source Catalog. These 
colors and magnitudes are listed in Paper VI. Although 
these galaxies likely have complex star formation histo- 
ries that complicate the determination of M/Lb, we can 
obtain a crude estimate by comparing their colors to sim- 
ple stellar population model grids. We use the Bruzual & 
Chariot (2003) models to obtain an average luminosity- 
weighted M I Lb for each galaxy. We then create seven 
logarithmic bins of mass and nonparametrically decom- 
pose the GC populations. These values are listed in Ta- 
ble □ 

Figure El shows the relationship between the metal- 
licity of the GC subpopulations and galaxy stellar mass. 
Similar to the previous figure, the metallicities of both 
GC subpopulations correlate strongly with galaxy mass. 
Across nearly four orders of magnitude in stellar mass, 
we find that Z oc Aff where (3 = 0.17 ± 0.04, 0.22 ± 



0.03, and 0.41 ± 0.01 for the metal-poor, metal-rich, and 
total GC populations. 

The more obvious manifestation of the correlation in 
the metal-poor GCs is due to the nature of the (g-z)- 
[Fe/H] relation, which steepens for bluer/metal-poor stel- 
lar populations. One implication of this is that it is a 
difficult task to tease out any differences between metal- 
poor clusters using broadband colors, and that the appar- 
ent universality of the population (at least in metallicity) 
may be in part due to photometric errors, small numbers 
of GCs, or the use of a less sensitive color such as V-I. 
Almost all previous studies on this topic have also used 
linear color-metallicity relations in either V-I or C-T\. 

Larsen et al. (2001) derive slopes between V-I and Mb 
of -0.016 ± 0.005 and -0.020 ± 0.008 for the blue and 
red GCs, respectively. If we convert these to [Fe/H] using 
the relation of Barmby et al. (2000), which has a slope of 
4.22, then this gives a — 0.17 and 0.21 for the blue and 
red GCs. As with the slopes we derived with KMM on 
individual galaxies, the blue GCs have a slightly steeper, 
and the red GCs a slightly shallower slope than what we 
finally derive with the binned sample. However, given 
the differences in sample, filters, and the assumption of 
a linear color-metallicity relation, the comparison can be 
considered consistent. 

We emphasize, however, that the significances of the 
correlations are independent of the color-metallicity re- 
lation. It is the exact values of these proportionalities, 
especially the one for the metal-poor GCs, that are crit- 
ically dependent on the adopted (<?-z)-[Fe/H] relation. 
We illustrate this in Figure El m which we plot the vari- 
ous linear fits for the Mg-[Fe/H] data assuming different 
transformations from (g-z) to [Fe/H]. The models we use 
are the same as those in Figure El an d as in that Fig- 
ure we linearly interpolate the models. While the slopes 
for the metal-rich GCs are all very similar, the slopes 
for the metal-poor GCs can vary by a factor of three. 
However, one bit of independent reassurance is visible in 
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Fig. 15. — Linear fits to the [Fe/H]-Afg relations assuming dif- 
ferent transformations from (g-z)o to [Fe/H]. We show fits using 
different models with different IMFs as well as our empirical broken 
linear fit to the Galactic GC data. While the slope of the relation 
in the red GC metallicity regime is generally robust across differ- 
ent assumed transformations, the slope of the blue GC relation can 
vary by a factor of three. 

Figures ITTA and ITBI where the large asterisk represents the 
Milky Way spheroid (bulge and halo), and its associated 
metal-rich and metal-poor populations of globular clus- 
ters. We take the luminosity of the Galactic spheroid to 
be M B = -19.19 from Cote (1999). That these values 
are close to what we might predict from the relations de- 
rived for Virgo ellipticals is interesting in itself, and lends 
credibility to our empirical color-metallicity transforma- 
tion. 

4.3. The Formation and Evolution of Globular Cluster 

Systems 

Our ability to characterize the GC metallicity distri- 
butions across a large range of galaxy masses provides a 
clearer and more precise view of the nature of globular 
cluster systems. The continuity of GC system proper- 
ties that we see paints a picture in which the formation 
and evolution of GC systems shares a common mecha- 
nism across all galaxies. Other properties of GC systems 
are also either constant or slowly varying, such as the 
luminosity function, the size distribution (Jordan et al. 
2005a), and the GC formation efficiency — globular clus- 
ters make up 0.25% of the baryonic mass of a galaxy 
(McLaughlin 1999). This continuity mirrors the struc- 
tural properties of the galaxies themselves, which show 
them to be part of a continuously varying single family 
(Ferrarese et al. 2005, Paper VI; Cote et al. 2005, Paper 
VIII). 

The correlations of the mean metallicities of metal- 
poor and metal-rich GCs with galaxy mass, luminos- 
ity, and color show that the formation of both sub- 
populations are closely linked to their parent galaxies. 
The nearly universal presence of a metal-poor popula- 
tion points to an era of star formation that was ubiq- 
uitous in the early universe at or shortly after reioniza- 
tion. Even at that early epoch — one associated with an 



epoch of "halo formation" — the metallicities of the form- 
ing star clusters are directly affected by the depth of the 
potential well that will eventually host the GC system. 
The metal-rich GCs appear to be associated with the for- 
mation of bulges and ellipticals — the metal-rich spheroid 
that defines the Hubble sequence in the local universe. 
Their numbers correlate strongly with the mass of the 
spheroid, indicating that they may have formed in the 
same events. Comparisons between GC and planetary 
nebula kinematics also support this view (Peng et al. 
2004). 

The relationship between the mean metallicity of a 
stellar system and its mass is a record of the enrichment 
and merging that occurred. There is a well established 
mass-metallicity relationship for the field stars of dwarf 
galaxies (e.g. Dekel & Woo 2003) and for the gas phase 
abundance of normal galaxies (Tremonti et al. 2004). Su- 
pernova wind feedback and its associated loss of metals 
has been used to explain the Z oc M aA relationship for 
Local Group dwarfs (Dekel & Silk 1986) and can repro- 
duce the flattening of this relationship at higher mass. 
In their high-resolution simulations, Kravstov & Gnedin 
(2005) find their simulated galaxies to have Z oc M 5 . 
Although we know little about the epoch of halo forma- 
tion, it is possible that galaxy outflows also play a role in 
enriching and triggering the formation of globular clus- 
ters (Scannapieco, Weisheit & Harlow 2004). Do these 
relationships apply to GC systems, and should they? 

Interestingly, we find Z oc M 0A1 for the total GC sys- 
tems of our sample galaxies. This may just be a coinci- 
dence since the slope is driven by the fraction of metal- 
rich GCs and hence must flatten at lower masses to fol- 
low the metal-poor relation. Moreover, it is in the dwarf 
regime where a = 0.4 is seen for the field stars, not for 
more massive galaxies. However, if we assume that the 
metal-poor and metal-rich GCs were formed in differ- 
ent events then we might expect that their individual a 
values would reflect the global metallicities of their host 
galaxies. In the mass range that we are concerned with 
(10 9 to 1O 12 M ), the slope of the M+-Z relationship for 
galaxies is significantly shallower than a — 0.4. Fitting 
the effective yield data for the SDSS galaxies presented 
in Table 4 and Figure 8 of Tremonti et al. (2004) , we find 
that in this mass range Z oc M 019 , which is quite simi- 
lar the relation for both GC subpopulations, and is also 
significantly shallower than the relation seen for Local 
Group dwarf galaxies. This suggests that star forma- 
tion and feedback in both the epochs of halo and metal- 
rich spheroid formation may not have been too different 
from each other or from that observed in the present 
day. In addition, we find that the Milky Way's total 
spheroid falls neatly on our relations derived for Virgo 
ellipticals. While there are some dependencies between 
the two samples because our color-metallicity relation 
was partly calibrated using Galactic GCs, the luminosity 
of the spheroid is independent. This agreement suggests 
that spheroids and GC systems in disk galaxies may form 
in much the same was as cluster ellipticals, although a 
large census of GC systems in disk galaxies will be nec- 
essary to explore such issues. 

The slopes of the mass-metallicity relationships are 
surprisingly similar for metal-rich and metal-poor GCs. 
One consequence of this is a nearly constant offset be- 
tween the metallicities of the two populations across 
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nearly three orders of magnitude in mass, A[Fe/H] oc 
M°- 05 ~ 1 dex. If the two populations are formed in 
different star forming events, this offset could point to a 
characteristic enrichment that is attained between star- 
bursts. 

Recent studies of galaxies in general point to the exis- 
tence of a bimodality in galaxy properties about a "char- 
acteristic mass" of 3 x 10 10 M Q (Kauffmann et al. 2003). 
Galaxies with higher masses tend to be spheroids with 
old stellar populations, and those with lower masses are 
likely to be star- forming disks. This bimodal rather 
than continuous distribution of galaxy properties may 
have its roots in the details of gas inflow and outflow, 
or feedback from active galactic nuclei (Dekel & Birn- 
boim 2005). For the GC systems, the only property that 
appears to change at this mass scale is the fraction of 
red GCs, which is nearly constant above this mass, and 
drops quickly below it. This reinforces the idea that the 
metal-rich GCs are linked to the formation of the metal- 
rich spheroid. This characteristic mass is also the regime 
at which the mass-metallicity relation of star forming 
galaxies appears to flatten (Tremonti et al. 2004). It is 
interesting to note then, although perhaps not entirely 
surprising since our sample consists only of early-types, 
that this characteristic mass does not appear to be im- 
portant for the mass-metallicity relationships of GC sys- 
tems. All galaxies appear to possess some old, GC-like 
stellar population, regardless of mass, gas content, or 
Hubble type (see Chandar, Whitmore, & Lee 2004, Olsen 
et al. 2004 for GCs in spiral galaxies and Seth et al. 2004 
for dwarf irregulars). More observations are necessary 
to determine whether the GC systems of disk galaxies 
other than the Milky Way are consistent with the same 
trends we see in this study, but if they are then it would 
point to a scenario where the processes that are affected 
by this characteristic mass scale occur after (or indepen- 
dently from) the formation of the GC system. This is 
further supported by the constant offset in metallicity 
between GCs and the field stars of ~ 0.8 dex that has 
been observed across a large range of galaxy luminosity 
(Jordan et al. 2004c). We speculate that this could be 
because massive star clusters form early in a star for- 
mation episode and are less affected by the subsequent 
feedback-related effects that shape the main body of the 
galaxy. 

With larger numbers of GCs in our sample allowing 
us to study galaxies in bulk, we choose to defer debates 
over whether a single GC system is unimodal, bimodal, 
or multimodal. Occasional claims for multimodality have 
been made for specific galaxies (in particular because of 
a population of intermediate-age metal-rich clusters), al- 
though it is very difficult to tease out such effects with 
single color data alone. In our data, we assess that VCC 
881 and 798 show color distributions that are potentially 
multimodal, but the occurrence is otherwise rare or too 
difficult to discern. While there are always individual 
galaxies with unique GC systems, the general trends are 
unmistakable. GC color distributions are on average bi- 
modal or asymmetric down to the magnitude limit of our 
sample, Mb ~ —15. Almost all galaxies have a popula- 
tion of metal-poor GCs, the mean metallicity of which 
increases with galaxy mass. These galaxies also possess 
some number of more metal-rich GCs whose number frac- 
tion and mean metallicity are also strong functions of 
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Fig. 16. — Metallicity dispersions (68% half-width) for the red 
and blue GC populations. Dispersions were calculated using the 
color dispersions and transforming them to metallicity using the 
proper branch of the (g— z)-[Fe/H] relation. Note that although 
the red GCs had a larger dispersion in color ( Figured, the steeper 
slope of the color-metallicity relation in the blue results in the blue 
GCs having the larger dispersion in metallicity. 



galaxy mass. 

The metallicity dispersions of GC systems increases 
rapidly for more massive galaxies, and this is largely 
driven by the increasing fraction of the red GCs. When 
comparing the metallicity dispersions of the red and blue 
GCs, the red GCs in massive galaxies do seem to have 
larger dispersions in color than the blue GCs, a factor 
of 1.3 in the mean. However, the shallower slope (i.e. 
increased metallicity sensitivity) of the (g-z)-[Fe/H] re- 
lation at [Fe/H] > —0.8 means that they do not translate 
into larger dispersions in metallicity. In fact, the slope 
of the color-metallicity relation is 2.8 ± 0.5 times steeper 
for metal-poor clusters, which means that it is the metal- 
poor GCs that have the larger dispersion in metallicity. 
This is illustrated in Figure where the median dis- 
persion is 0.53 ± 0.07 dex for the metal-poor GCs, and 
~ 0.29±0.04 dex for the metal-rich GCs, where the errors 
reflect the formal errors in the slope of the (g-z)-[Fe/H] 
relation. 

The dispersion in metallicity potentially tells us some- 
thing about the degree of gas mixing that occurred 
during galaxy formation. A large metallicity disper- 
sion could be caused by prolonged formation in a single 
galaxy, or also by a combination of star formation events 
that are chemically isolated from each other (whether 
spatially or temporally) and that only merge together af- 
ter star formation has ceased. This latter scenario could 
apply to both metal-poor globular clusters that formed 
from isolated gaseous fragments in the early universe, or 
to metal-rich clusters that form during the hierarchical 
build up of the metal-rich spheroid. It will be interesting 
to see whether the lack of a larger dispersion in the metal- 
rich GCs, if confirmed, poses a problem for secondary 
formation scenarios that require the red GCs to be the 
product of multiple star forming events at moderate to 
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low redshift, or if the metallicity resolution is still inad- 
equate to resolve these effects. A narrower dispersion in 
metallicity also suggests that the metal-rich population 
may be formed in few large starforming events where the 
gas is well mixed. 

Naturally, this inference is highly dependent on the 
color-metallicity relation, but we point out that none of 
the models shown in Figure IT21 will produce a larger dis- 
persion for the metal-rich GCs. At best, the two pop- 
ulations have the same dispersion. The Milky Way GC 
system does not show such a large dispersion difference 
between the metal-poor and metal-rich GCs, however, 
and so we will likely need the final (g-z)- [Fe/H] calibra- 
tion before we can verify this intriguing result. 

5. CONCLUSIONS 

Ever since the launch of HST, the color distributions 
of extragalactic globular cluster systems have provided a 
unique window onto the evolution of elliptical galaxies. 
In this paper we present the globular cluster color dis- 
tributions in (g-z) of 100 early- type galaxies from the 
ACS Virgo Cluster Survey. As the largest sample of 
its kind, studied with the increased sensitivity and spa- 
tial sampling of the ACS/WFC, and with the 3475 and 
Z850 bandpasses, this is one of the definitive data sets for 
extragalactic GC studies. We use carefully constructed 
control fields to account for galaxy specific and spatially 
varying foreground and background contamination. We 
also use a new size-magnitude algorithm to produce a 
clean sample of GCs. Using both two-Gaussian models 
and a nonparametric decomposition of the color distri- 
butions, we measure the characteristics of the traditional 
metal-poor and metal-rich subcomponents of the GC sys- 
tems. We also present a preliminary transformation from 
(g-z) to [Fe/H] using the colors of Galactic GCs, as well 
as GCs in M49 and M87. Our main conclusions are: 

1. While the color distributions of individual galaxies 
can show significant variations from one another, 
their general properties are consistent with con- 
tinuous trends across galaxy luminosity, color, and 
mass. Both dwarf and giant ellipticals appear to be 
subject to the same family of processes that drive 
the formation and evolution of GC systems. 

2. Galaxies at all luminosities in our study, on av- 
erage, appear to have bimodal or asymmetric GC 
color distributions. All galaxies have a system of 
blue GCs with an average color around g-z ~ 0.9 
and a red component whose fraction varies with the 
luminosity or color of the galaxy. When decom- 
posed into blue and red GC subpopulations, the 
red GCs on average compose ~ 15% of the GCs 
in the faintest and bluest galaxies in our sample 
(M B = -16, g-z = 1.0), and up to 60% of the GCs 
the brightest and reddest galaxies (Mb — —21, g- 
z = 1.5). 

3. The colors of the blue and red GC subpopulations 
correlate with the luminosity and color of the host 
galaxy. In (g-z), the slope determined for the red 
GCs is 4.6 times steeper than for the blue GCs. 

4. The color and metallicity widths of the entire GC 
system increases with galaxy luminosity and color. 



The width of the red population is larger than that 
of the blue population in color, but because of the 
decreased sensitivity of color to metallicity in the 
blue, the widths of the red peak in metallicity are 
consistent with being smaller than those of the blue 
peak. 

5. Using new Galactic GC imaging combined with 
M49 and M87 GCs that have spectroscopic metal- 
licities, we have defined a new preliminary color- 
metallicity relation for globular clusters. We find 
that this relation shows clear nonlinearity and 
parametrize the relation with a broken line with 
a break point at (g-z)=1.05, [Fe/H]= -0.81. The 
steepening of the relation at low metallicity is crit- 
ical to the proper interpretation of metal-poor GC 
colors. Current evolutionary synthesis models pro- 
duce colors for 13 Gyr SSPs that are roughly con- 
sistent with our color-metallicity relation, but show 
large variation for g-z < 1.0. 

6. The metallicities of the metal-rich and metal-poor 
GC subpopulations both correlate with the lumi- 
nosity of the host galaxy, with the slope for the 
metal-rich GCs 1.7 times steeper. The derived 
slopes are not nearly as disparate in metallicity 
as they are in color because of the steeper color- 
metallicity relation for metal-poor GCs. How- 
ever, uncertainty in the (g— z)— [Fe/H] relation in 
the blue can cause this slope to change by a factor 
of three. We find that Z cx L 016±0 04 for metal- 
poor GCs, Z oc LO.25io.03 for m etal-rich GCs , and 
Z cx iO.52io.02 for the tQtal GC p p U i a tion. 

7. We use the optical-infrared colors of the ACSVCS 
galaxies to derive crude stellar masses and find 
that the metallicities of the GC populations also 
correlate with the stellar mass of the host galaxy. 
The ratio of the slopes is smaller again, with the 
slope for the metal-rich GCs 1.4 times steeper. We 
find Z oc M° 17±0 04 for metal-poor GCs, Z oc 
M o.22io.03 for me tal-rich GCs, and Z oc M° A1±om 
for the total GC population. 

8. The small difference in slope in the M+-Z relation 
means that the metallicity difference between the 
metal-poor and metal-rich GC subpopulations is a 
nearly constant A [Fe/H] ~ 1 dex, varying slowly 
with stellar mass as M 05 . 

As part of the ACSVCS, we have also obtained long-slit 
spectroscopy of the entire galaxy sample. In a subsequent 
paper, we will present the properties of the GC color 
distributions as a function of galaxy metallicity, alpha- 
enhancement, and total dynamical mass. A future paper 
will also address these results in the context of detailed 
modeling of host galaxy merging histories. 
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TABLE 3 

Linear fits to Red and Blue GC components of individual galaxies versus M b 





p- value 


(g - z) = a + b x M B 
a 


b 


Red GCs 


p < 0.05 


0.629 ±0.082 


-0.036 ± 0.004 




< p < 1 


0.503 ±0.061 


-0.042 ± 0.003 


Blue GCs 


p < 0.05 


0.424 ± 0.063 


-0.026 ±0.003 




< p < 1 


0.462 ± 0.046 


-0.024 ±0.002 



TABLE 4 

Two-Component Decompositions for GC Color Distributions, Binned by M b 



M b 




fir 


s 68,b 


«68,r 

(5) 


fred 


((g-z)host) 


Ngc 


Mi 




(1) 


(2) 


(3) 


(4) 


(6) 


(7) 


(8) 


(9) 


(10) 


-21.4 


0.95 ± 0.01 


1.40 ±0.01 


0.13 ±0.01 


0.17 ±0.01 


0.63 ±0.03 


1.50 


3214 


1.26 ±0.01 


0.30 ±0.01 


-20.5 


0.93 ± 0.02 


1.32 ±0.02 


0.12 ±0.01 


0.16 ±0.01 


0.56 ± 0.06 


1.47 


3122 


1.15 ±0.01 


0.26 ±0.01 


-19.5 


0.91 ± 0.01 


1.28 ±0.02 


0.10 ±0.01 


0.16 ±0.01 


0.52 ± 0.04 


1.44 


1692 


1.08 ±0.01 


0.26 ±0.01 


-18.5 


0.89 ± 0.02 


1.24 ±0.03 


0.10 ±0.01 


0.16 ±0.01 


0.51 ±0.06 


1.38 


1088 


1.04 ±0.01 


0.24 ± 0.01 


-17.7 


0.89 ±0.01 


1.16 ±0.02 


0.10 ±0.01 


0.12 ±0.01 


0.33 ± 0.04 


1.36 


644 


0.96 ±0.01 


0.18 ±0.01 


-16.6 


0.87 ±0.01 


1.13 ±0.03 


0.10 ±0.01 


0.12 ±0.01 


0.29 ±0.06 


1.14 


505 


0.92 ±0.01 


0.17 ±0.01 


-15.7 


0.88 ± 0.01 


1.12 ±0.04 


0.10 ±0.01 


0.09 ± 0.02 


0.15 ±0.06 


1.14 


233 


0.89 ±0.01 


0.14 ±0.01 



Mean galaxy magnitude in bin 

Median colors of blue and red GC components 
68% half- widths of blue and red GC components 
Fraction of clusters in red component 
Mean color of galaxies in this bin 
Number of GCs 

Median of entire GC distribution 

68% half-width of entire GC distribution 



TABLE 5 

Two-Component Decompositions for GC Color Distributions, Binned by galaxy (g - z) 



(9-2)0 

(1) 


(2) 


(3) 


•568,6 
(4) 


S68,r 

(5) 


fred 

(6) 


Ngc 

(7) 


Mi 

(8) 


(9) 


1.52 


0.95 ±0.02 


1.40 ±0.02 


0.13 ±0.01 


0.16 ±0.01 


0.55 ± 0.05 


2064 


1.20 ±0.01 


0.30 ±0.01 


1.47 


0.91 ±0.00 


1.32 ±0.01 


0.10 ±0.01 


0.18 ±0.01 


0.58 ± 0.02 


6375 


1.13 ±0.01 


0.28 ±0.01 


1.37 


0.91 ±0.02 


1.23 ±0.02 


0.12 ±0.01 


0.15 ±0.01 


0.47 ±0.05 


1439 


1.04 ±0.01 


0.23 ±0.01 


1.26 


0.88 ± 0.02 


1.16 ±0.05 


0.10 ±0.01 


0.13 ±0.02 


0.31 ±0.08 


518 


0.94 ± 0.01 


0.18 ±0.01 


1.14 


0.89 ±0.01 


1.12 ±0.03 


0.11 ±0.01 


0.09 ±0.01 


0.18 ±0.05 


393 


0.92 ±0.01 


0.14 ±0.01 


1.02 


0.87 ±0.01 


1.14 ±0.06 


0.13 ±0.02 


0.10 ±0.04 


0.14 ±0.07 


99 


0.88 ±0.02 


0.15 ±0.02 



Mean galaxy color in bin 

Median colors of blue and red GC components 
68% half-widths of blue and red GC components 
Fraction of clusters in red component 
Number of GCs 

Median of entire GC distribution 

68% half-width of entire GC distribution 
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TABLE 6 

Linear fits to Red and Blue GC components of versus Galaxy Properties 



X 

Galaxy Property 


Y 

GCS Property 


GCs 


Y = a + bX 
a 


b 




M B 


{g-z) 


Red 


0.163 ±0. 


006 


-0.058 ± 0. 


.001 






Blue 


0.668 ± 


.046 


-0.013 ±0 


.003 


(g-z) 


{g-z) 


Red 


0.405 ± 0. 


091 


0.622 ± 0. 


.063 


Blue 


0.795 ± 0. 


033 


0.078 ±0. 


024 


M B 


[Fe/H] 


Red 


-2.452 ±0. 


021 


-0.105 ±0 


.013 




Blue 


-2.771 ± 0. 


232 


-0.064 ± 0. 


.015 






All 


-4.854 ±0. 


124 


-0.206 ± 0. 


006 


loguiiM*) 


[Fe/H] 


Red 


-2.759 ± 0. 


.032 


0.219 ±0. 


028 






Blue 


-3.296 ±0. 


328 


0.167 ±0. 


039 






All 


-5.250 ± 0. 


156 


0.409 ± 0. 


014 



TABLE 7 

Two-Component Decompositions for GC Color Distributions, Binned by Stellar Mass 



M*(M ) 
(1) 


(2) 


fir 

(3) 


s 68,f> 
(4) 


•568, r 

(5) 


fred 

(6) 


{{g-z)ho S t) 

(7) 


Ngc 

(8) 


Ml 
(9) 


(10) 


4.8E+11 


0.96 ±0.02 


1.40 ±0.01 


0.14 ±0.01 


0.16 ±0.01 


0.58 ± 0.04 


1.50 


3558 


1.23 ± 0.01 


0.30 ± 0.01 


1.7E+11 


0.94 ±0.02 


1.32 ± 0.01 


0.12 ±0.01 


0.16 ±0.01 


0.57 ±0.04 


1.48 


2789 


1.16 ± 0.01 


0.26 ±0.01 


5.7E+10 


0.89 ±0.01 


1.26 ±0.02 


0.10 ±0.01 


0.16 ±0.01 


0.54 ± 0.04 


1.46 


1326 


1.07 ±0.01 


0.25 ±0.01 


2.0E+10 


0.92 ±0.02 


1.27 ±0.02 


0.12 ±0.01 


0.15 ±0.01 


0.44 ± 0.05 


1.39 


1061 


1.04 ±0.01 


0.24 ±0.01 


8.6E+09 


0.88 ±0.01 


1.16 ±0.02 


0.10 ±0.01 


0.13 ±0.01 


0.37 ±0.04 


1.33 


528 


0.96 ±0.01 


0.18 ±0.01 


2.7E+09 


0.88 ±0.02 


1.13 ±0.04 


0.10 ±0.01 


0.12 ±0.01 


0.32 ±0.07 


1.20 


487 


0.93 ±0.01 


0.16 ±0.01 


1.1E+09 


0.87 ±0.01 


1.13 ±0.04 


0.11 ±0.01 


0.09 ±0.03 


0.16 ±0.06 


1.11 


268 


0.88 ±0.01 


0.15 ±0.01 



''--n galaxy stclla 



2,3, 
5 

68% half-widths of blue and red GC components 



Median colors of blue and red GC components 

6. 



Fraction of clust 

7 



Mean color of galaxies in this bin 
Number of GCs 

Median of entire GC distribution 

68% half-width of entire GC distribution 



